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ABSTRACT 

We investigate, using the modeling code SHAPE, the three-dimensional 
structure of the bipolar Homunculus nebula surrounding Eta Carinae as 
mapped by new ESO VLT/X-Shooter observations of the H 2 A = 2.12125)Rm 
emission line. Our results reveal for the first time important deviations 
from the axisymmetric bipolar morphology: 1) circumpolar trenches in 
each lobe positioned point-symmetrically from the center and 2) offplanar 
protrusions in the equatorial region from each lobe at longitudinal 55°) 
and latitudinal ( 10 ° — 20 °) distances from the projected apastron direction 
of the binary orbit. The angular distance between the protrusions (~ 110°) 
is similar to the angular extent of each polar trench 130°) and nearly 
equal to the opening angle of the wind- wind collision cavity 110°). As 
in previous studies, we confirm a hole near the centre of each polar lobe 
and no detectable near-IR H 2 emission from the thin optical skirt seen 
prominently in visible imagery. We conclude that the interaction between 
the outflows and/or radiation from the central binary stars and their ori- 
entation in space has had, and possibly still has, a strong influence on the 
Homunculus. This implies that prevailing theoretical models of the Ho- 
munculus are incomplete as most assume a single star origin that produces 
an axisymmetric nebula. We discuss how the newly found features might be 
related to the Homunculus ejection, the central binary and the interacting 
stellar winds. 

Key words: methods: numerical - stars: individual: Eta Carinae - circumstellar 

matter - stars: mass-loss - stars: winds, outflows 


1 INTRODUCTION 


* Based on observations made with ESO Telescopes at 
the La Silla Paranal Observatory under programme ID 
088.D-0873(A) (PI Groh). 
t e-mail: wsteffen@astro.unain.mx 
t CNPq/Science without Borders Fellow 
§ NASA Postdoctoral Program Fellow 


Massive stars have major impacts on their host galaxies via 
the input of ionizing photons, energy, and momentum into 
the interstellar medium. They also contribute to the chem- 
ical evolution of the Universe via the ejection of enriched 
elements through stellar winds and supernovae. Our cur- 
rent understanding is that some massive stars undergo vio- 
lent ejections of a fraction to several-tens of solar masses as 
they evolve through the luminous blue variable (LBV) stage 
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(Langer et al. 1994; Humphreys & Davidson 1994; Smith & 
Owocki 2006). The mechanisms responsible for these erup- 
tive events, as well as their frequency and the total amount 
of mass lost, are not currently understood. 

Eta Carinae is the closest and most studied example of 
an object that underwent a ‘Great Eruption’ (Davidson & 
Humphreys 2012). At least lOM©, and potentially more than 
4OM0, were ejected during this event (Smith et al. 2003a; 
Smith 2006; Gomez et al. 2010), causing the formation of 
a dusty bipolar nebula known as the Homunculus (Gaviola 
1950). Analysis of the proper motion of the Homunculus 
shows that the eruption happened in the 1840s (Gurrie et 
al. 1996; Morse et al. 2001). This, combined with its relative 
proximity (2.3 kpc, Allen & Hillier 1993; Smith 2006), means 
that the Homunculus’s structure can be investigated in an 
unparalleled level of detail. This provides valuable insights 
into the mechanisms responsible for LBV eruptions and their 
role in the evolution of the most massive stars. 

Eta Garinae itself is a massive, highly eccentric (e ~ 0.9) 
binary system (Damineli et al. 1997) with a combined lu- 
minosity and mass of L* ~ 5 x 10® Lq (Gox et al. 1995; 
Davidson & Humphreys 1997) and M* > 110 M© (Hillier 
et al. 2001), respectively. The LBV primary is responsible 
for most of the system luminosity (and thus mass), with 
a mass-loss rate of 8.5 x Mq yr~^ and wind termi- 
nal velocity of ~ 420 kms“^ (Groh et al. 2012a; see also 
Hillier et al. 2001, 2006). This extremely dense wind causes 
the photosphere to be formed in the wind, and one infers 
an effective temperature of Tqe ~ 9400 K at optical depth 
T = 2j3 (Hillier et al. 2001, 2006; Groh et al. 2012a). Only 
indirect constraints are available for the companion star, 
which likely has a temperature of Tefr ~ 36, 000 — 41, 000 K, 
luminosity 10® L© ^ T* < 10® L© (Mehner et al. 2010), 
wind terminal speed Voo ~ 3000 kms“^, and mass-loss rate 
M ~ 1.4 X 10“® Mq yr“^ (Pittard & Gorcoran 2002; Parkin 
et al. 2009). Gurrently, the companion has a significant im- 
pact on the wind of the primary (Pittard & Gorcoran 2002; 
Okazaki et al. 2008; Parkin et al. 2009; Madura & Groh 2012; 
Groh et al. 2010a,b, 2012a,b; Madura et al. 2012, 2013). 
However, it is unclear whether the companion had any role 
in triggering the eruption or shaping the Homunculus. 

Numerous theoretical studies have investigated the for- 
mation of the Homunculus. These studies fall mainly into 
six categories: interacting winds (Prank et al. 1995, 1998; 
Dwarkadas & Balick 1998; Langer et al. 1999; Gonzalez et 
al. 2004a,b), binary interaction (Soker 2001, 2004, 2007; 
Kashi & Soker 2010), a radiatively-driven wind from a 
rapidly-rotating star (Owocki & Gayley 1997; Maeder & 
Desjacques 2001; Dwarkadas & Owocki 2002; Owocki 2003, 
2005), a thermally-driven magnetohydrodynamic (MHD) ro- 
tating wind (Matt & Balick 2004), explosive mass loss from 
a rapidly-rotating star (e.g. Smith & Townsend 2007), and 
stellar mergers (Gallagher 1989; Iben 1999; Morris & Pod- 
siadlowski 2006). The interacting winds scenario requires a 
pre-existing slow and dense equatorial torus to pinch the 
waist of the bipolar structure. However, observations show 
that the thin equatorial disc or ‘skirt’ seen in optical im- 
ages of the Homunculus is no older than the bipolar lobes 
and that there is not enough mass in the skirt to shape 
the overall morphology of the nebula (Smith & Gehrz 1998; 
Davidson et al. 2001; Morse et al. 2001; Smith 2006). The 
MHD model of Matt & Balick (2004) requires a huge field 


strength (> 2.5 x 10^ G at the equator) to shape the large 
amount of mass ejected, and it is unclear if such a field was 
available. Models employing radiatively-driven winds from 
a rapidly-rotating star have difficulties reproducing simulta- 
neously the bipolar structure and the thin equatorial skirt. 
The large amount of mass expelled over such a short period 
of time also cannot be readily explained by the standard line- 
driven-wind formalism, and may even surpass the capabil- 
ity of a super-Eddington continuum- driven wind (Owocki et 
al. 2004; Smith & Townsend 2007). However, a continuum- 
driven explosion from a rapidly-rotating single star may ex- 
plain both the bipolar shape and thin equatorial skirt (Smith 
& Townsend 2007). Stellar mergers have difficulty account- 
ing for previous and subsequent large eruptions, such as the 
smaller eruption in 1890 that formed the Little Homunculus 
(Ishibashi et al. 2003; Smith et al. 2004). While mergers may 
seem unlikely, they also cannot yet be fully ruled out. 

Besides the single-star axisymmetric model of Smith & 
Townsend (2007), another possibility for explaining the for- 
mation and shape of the Homunculus involves some sort 
of binary interaction, perhaps during periastron passage 
(Kashi & Soker 2010; Smith 2011). Possible evidence for such 
a scenario is the correlation of periastron passage with the 
peaks observed in Eta Garinae’s light curve in the decades 
leading up to the Great Eruption (see fig. 2 of Damineli 
1996 and also Smith & Erew 2011). Moreover, it was re- 
cently found that the orbital axis (orthogonal to the or- 
bital plane and through the system centre of mass) of the 
Eta Garinae binary is closely aligned in three-dimensional 
(3D) space with the inferred polar symmetry axis of the 
Homunculus (Madura et al. 2012). Unfortunately, all cur- 
rent models for the Homunculus that involve binary inter- 
action are phenomenological and rely on assumptions such 
as strong accretion by the secondary star. These simple mod- 
els are incapable of predicting the exact shape and degree of 
asymmetry in any resulting nebula (Soker 2001). Numerical- 
hydrodynamical simulations are required. 

If the binary companion had any influence on the for- 
mation or shaping of the Homunculus, it should have left 
an imprint in the Homunculus’s morphology. There are at 
least two ways the binary companion and Homunculus’s for- 
mation might be related, and both could be relevant. Eirst, 
the companion may be directly responsible for triggering 
the eruption itself, say via a stellar collision (Smith 2011). 
Second, the companion could influence the shape of the neb- 
ula during and/or after any explosion, via a collimated fast 
wind like that proposed by Soker (2001), or simply via its 
extremely fast, low-density wind and orbital motion. 

Of course, any theoretical model for the Homunculus 
can only be as good as the observations on which it is based. 
Several observational studies, using mainly long-slit optical 
spectroscopy and mid-infrared imaging, have investigated 
the 3D structure of the Homunculus. Ground-based optical 
spectroscopy at medium spectral resolution has long been 
employed to measure doppler-shifts of emission lines formed 
in the Homunculus lobes, showing that the nebula is bipo- 
lar with an approaching south-east (SE) lobe, a receding 
north-west (NW) lobe, and a thin equatorial skirt (Thack- 
eray 1951, 1956a,b, 1961; Meaburn et al. 1987, 1993; Hillier 
& Allen 1992; Allen & Hillier 1991, 1993; Hillier 1997; Gur- 
rie et al. 1996; Meaburn 1999). The improved spatial res- 
olution of Hubble Spaee Teleseope (HaST) /S pace Telescope 
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Imaging Spectrograph (STIS) optical observations allowed 
further detailed studies of the shape and orientation of the 
Homunculus lobes (Davidson et ah 2001). The study of near- 
IR (NIR) emission lines discovered the presence of molecu- 
lar gas (H2; Smith & Davidson 2001), making it possible 
to trace the structure of the optically-obscured back-side of 
the Homunculus (Smith 2002, 2006; Teodoro et ah 2008). 
Smith et ah (2003a) proposed that the H2 gas is likely asso- 
ciated with the outer, cool dust shell seen in mid-IR images 
(Hackwell et ah 1986; Smith et ah 1995; Smith, Gehrz & 
Krautter 1998; Polomski et ah 1999; Chesneau et ah 2005), 
which traces the majority (~ 90%) of the Homunculus mass 
(15-45 M©; Smith et ah 2003a; Smith 2006; Gomez et ah 
2010). NIR observations (Smith 2002) also showed that the 
H2 V = 1 — 0 S(l) emission originates in a shell located ex- 
terior to the region of [Fe ii] and [Ni ii] emission investigated 
by Davidson et ah (2001). 


Subsequent ground-based, high-spectral-resolution ob- 
servations of H2 and [Fe ii] lines provided the most detailed 
view of the Homunculus to date (Smith 2006). This study 
confirmed the double-shell structure of the Homunculus, 
with inner [Fell] and outer H2 shells. It also showed that 
about 75% of the Homunculus mass is located at high lat- 
itudes (above 45°). In addition, more than 90% of the ki- 
netic energy was released at these high latitudes during the 
Great Eruption, which seems to have had a short duration 
(< 5 years) given the thinness of the H2 shell (Smith 2006). 
However, these observations had limited spatial coverage, 
using only five long-slits aligned along the major polar axis 
of the Homunculus. This provides insufficient information 
for a complete 3D study of the structure of the Homuncu- 
lus. Moreover, all previous observational studies of the Ho- 
munculus focused on its large-scale axisymmetric shape in 
the context of a single-star model. Departures from axisym- 
metry and the role of the companion are either only briefly 
mentioned or neglected. Detailed knowledge of the 3D struc- 
ture of the Homunculus, including small-scale departures 
from axisymmetry, is crucial for any proper theoretical mod- 
eling of its formation and the Great Eruption. 


Using an extensive new set of ESO Very Large Telescope 
(VLT)/X-shooter spectral mapping observations of the en- 
tire Homunculus, we investigate the Homunculus’s full 3D 
structure, providing a more detailed view of its shape. In 
particular, we probe deviations from axisymmetry to deter- 
mine whether the companion star might have played a role 
in shaping the Homunculus. We analyze the X-shooter data 
cubes using the SHAPE software (Steffen et al. 2011) and 
generate, for the first time, a complete 3D model of the Ho- 
munculus that includes smaller-scale features. In this initial 
investigation we focus on the H2 v = 1 — 0 S(l) emission 
line, which traces the fronts and backs of the Homunculus’s 
polar lobes (Smith 2006). 


This paper is organized as follows. In Sect. 2 we present 
our new observational mapping data obtained with VLT /X- 
shooter. Sect. 3 describes the SHAPE modeling procedure 
and Sect. 4 the results. We discuss the results in Sect. 5 
before summarizing our conclusions in Sect. 6. 



5 0-5 
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Figure 1. Mapping position of each slit superimposed over an 
HST/ ACS F550M image of the Homunculus obtained on 4 Au- 
gust 2006 (HST Proposal ID 10844, PI: K. Davidson; Martin et 
al. 2006). The original image has been rotated by 49° counter- 
clockwise and corrected for the mean expansion of the nebula 
between 2006.6 and 2012.2 0.20" (Smith & Gehrz 1998). Hor- 

izontal red lines mark the numbered slit locations where select 
position-velocity (P-V) images used in this paper were extracted 
from the observational data cube (see Fig. 3). 


2 OBSERVATIONS, DATA REDUCTION AND 

DATA CUBE ASSEMBLING 

A full spectral mapping (2990 A through 24790 A) of the 
entire Homunculus nebula (14" x 20") was obtained on the 
nights of 9 and 10 March 2012, using the ESO VLT/X- 
shooter spectrograph (ESO program ID 088.D-0873(A)). A 
detailed description of the instrument’s capabilities and con- 
figurations is given in Vernet et al. (2011). Since we focus 
on the H2 u = 1 — 0 S(l) A = 2.12125 yum^ emission line to 
model the Homunculus, we discuss here only the parameters, 
processing, and reductions for the X-shooter NIR arm. 

The data were obtained using the 0.4" x 11" slit, which 
resulted in a spectrum with R = X/6X = 11300 and a sam- 
pling of 2 pixels per EWHM. The mapping of the nebula was 
done with the slit oriented at position angle P.A.=— 41°, i.e. 
perpendicular to the projected major axis of the Homuncu- 
lus. A total of 92 dithered positions along the nebula was 
used in the full mapping (Eig. 1). The spatial coverage was 
from —6.93" to 4-6.93" from the central source in the X di- 
rection, and from —10.02" to 4-9.93" in the Y direction. The 
individual exposure times ranged from 0.06651 seconds close 
to the central source to 30 seconds in the lobes in order to 
not overexpose. Multiple images at each position were ob- 


^ All wavelengths in this paper refer to measurements in air. 
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tained and combined. The total exposure time thus ranged 
from 30 seconds close to the star to 150 seconds in the lobes. 

Data processing and reduction were performed using 
the X-shooter pipeline (cpl-6.1.1 and xsh/ 1.5.0; Goldoni 
et al. 2006) for which the standard processing includes re- 
moval of dark current, order identification and tracing, cor- 
rection for pixel-to-pixel sensitivity variations, dispersion so- 
lution, master response, and merging of all orders into a 
single wavelength calibrated 2D spectrum (for each slit po- 
sition). In the spectral direction, each reduced 2D spectrum 
ranges from 9940 A to 24790 A with a dispersion of 1 A per 
pixel, whereas in the spatial direction, the final scale for the 
NIR data is 0.21" per pixel. 

The assembling of the final data cube was achieved by 
creating a grid with spatial dimensions set by the maximum 
displacement between the offset positions along slit width 
and length. Each spatial pixel in the data cube grid was 
set to a size of 0.21" x 0.21". The size of each pixel in the 
spectral direction was not changed. Therefore, in the final 
data cube, each voxel has 0.21" x 0.21" x 1.0 A. In the case 
where slits overlap each other, we adopted the average value 
for each overlapping pixel. 

As an illustration of the high quality of the final assem- 
bled data cube. Fig. 2 shows H 2 iso- velocity maps for which 
the spectra were normalized by the adjacent continuum to 
highlight the structures in emission (corresponding to mate- 
rial in the Homunculus lobes). Although line emission from 
the backs of the lobes is weaker than that coming from the 
fronts (see Figs. 2e, /, and /), we were still able to map them. 


3 3D STRUCTURE MODEL 

With the 3D morpho-kinematic modeling we wish to chiefiy 
address two questions: 1) Does the small-scale structure of 
the Homunculus, beyond the cylindrical- symmetric compo- 
nent of earlier models, include information on the details of 
the ejection mechanism? 2) Are there symmetries that might 
indicate a relation to the central binary interaction? 

Previously constructed morpho-kinematic models of the 
Homunculus assumed axisymmetry, yielding only its overall 
shape (e.g. Davidson et al. 2001; Smith 2006; Teodoro et al. 

2008), although departures from true axial symmetry are 
known (Morse et al. 1998; Smith 2006). In this paper we con- 
sider three refinements to previous models. The first is the 
discrepancy between models based on Doppler- velocity mea- 
surements and those based on images. Second, we attempt 
to reconstruct some of the small-scale structures in the polar 
regions. Third, we include two newly- found kinematic fea- 
tures near the equatorial region of the Homunculus. As in 
Smith (2006), we find no detectable NIR H 2 emission from 
the thin equatorial skirt that is prominent in optical im- 
ages of the Homunculus. We therefore model only the main 
bipolar structure. 

We first consider the overall structure of the Homuncu- 
lus. Some papers present image- and velocity-based mod- 
els in the same work (e.g. Davidson et al. 2001). The sys- 
tematic discrepancy between models based on images and 
spectroscopy is that the spectroscopic reconstructions from 
molecular hydrogen yield a slimmer Homunculus than the 
image-based models. Several possible explanations come to 
mind for such a method-dependent difference in the recon- 


structions. First, the differences could be intrinsic, since the 
image-based reconstructions are based on continuum scat- 
tered by dust, whereas the Doppler-reconstructions use line 
emission that might originate in spatially distinct regions. 
Furthermore, the fact that Eta Carinae is a binary with a 
highly elliptical orbit and interacting winds could produce 
a symmetry that is non-cylindrical. 

The impression from the optical images, defined largely 
by dust-scattered continuum light from Eta Carinae interior 
to the Homunculus shell, is that a relatively thin, clumpy 
surface defines the location of the scattering dust. One is left 
with the impression that the bulk of the H 2 emission is just 
outside this surface. Molecular hydrogen is thus expected 
to be located mainly outside the dust shell and hence also 
outside the continuum emission. However, reconstructions 
from molecular line emission show a slimmer Homunculus 
than expected and hint toward a problem in the assumptions 
of the velocity field. 

Reconstructions based on Doppler-velocity measure- 
ments of expanding nebulae often assume homologous ex- 
pansion, i.e. radial expansion with the magnitude of the ve- 
locity increasing linearly with distance. This is a very strin- 
gent assumption and can be met only under certain con- 
ditions, and only to some level of accuracy. Hydrodynam- 
ical simulations of wind-blown bipolar nebulae show that 
substantial distortion can occur in 3D reconstructions if de- 
viations from a homologous expansion are not taken into 
account (Steffen, Garcia-Segura & Koning 2009). We there- 
fore explore whether hydrodynamically-plausible deviations 
from a homologous expansion can lead to a more consistent 
3D reconstruction of the large-scale structure. 

In addition to the large-scale structure, we attempt to 
reconstruct the most salient small-scale features based on 
the molecular hydrogen Doppler-velocity measurements pre- 
sented in Figs 2 and 3. Here it is important to note that we 
reconstruct a single continuous surface of each lobe and we 
assume that the deviations from a homologous expansion 
are only large-scale, such that the local structures are rea- 
sonably well mapped by the velocity measurements. 

As will be discussed later, some smaller-scale structures 
in the polar regions of the Homunculus may be caused by 
hydrodynamical instabilities, that at some time in the past 
must have introduced local velocity disturbances. If instabil- 
ities have been active in the relatively distant past and are 
now only passive, then the local deviations from homologous 
expansion will be small. Detailed numerical hydrodynamic 
simulations are needed to test the accuracy of the assump- 
tion of locally homologous expansion in the complex polar 
regions of the Homunculus. 

Assuming that deviations from a homologous expansion 
are only large scale, we first investigate the large-scale Ho- 
munculus shape and then the small-scale structure of the 
polar and equatorial regions. Afterward, we investigate the 
effects of dust in order to understand better the observed 
brightness distribution in the H 2 data cube and its repre- 
sentation as position-velocity (P-V) diagrams. 

3.1 Modeling method 

We constructed the 3D model using the morpho-kinematic 
and radiative transfer code SHAPE (Steffen et al. 2011). 
Based on images and spectroscopic data, SHAPE has been 
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Figure 2. Iso-velocity images of H 2 t’ = 1 — 0 S(l) at 2.12125 The structures seen in emission in each slice were mapped both in 
the spatial and spectral directions, allowing the reconstruction of the 3D shape of the emitting volume. The black horizontal lines across 
the images are regions where the spacing between adjacent slits was larger than the slit width and, therefore, contain no data. 


used extensively for modeling the complex 3D structure of 
planetary nebulae (e.g. Garcfa-Dfaz et al. 2012; Clark et 
al. 2013), novae (Ribeiro et al. 2013), radiation transfer in 
dusty pre-planetary nebulae (Koning et al. 2013), quark- 
nova ejecta (Ouyed et al. 2012), and the complex planetary 
nebula Hubble 5 (Lopez et al. 2012), which in shape is rather 
similar to the Homunculus. The model is constructed inter- 
actively with 3D mesh structures (see Fig. 4). A velocity field 
and a brightness distribution is then assigned. We use the 
physical renderer of SHAPE which allows radiative transfer 
with dust opacity. Radiative transfer is helpful for modeling 
of the small-scale structures in the Homunculus since the ef- 
fects of continuum absorption by dust cannot be neglected. 
Therefore, we introduce a dust shell that is basically a copy 
of the H 2 emission shell, but with a slightly smaller size so 
that they only partially overlap. 

For this work the details of the dust properties are not 
relevant, but for completeness we mention that we use op- 
tical properties appropriate for silicate dust with a power- 
law size distribution between 0.5 and 2.5 microns, a power- 
law index of 3.5, and optical properties according to Krugel 
(2008). The H 2 emissivity and dust density are assumed 
to be a smooth distribution that decreases slowly with dis- 


tance. Therefore, the resultant, often strong, projected im- 
age model brightness variations are mainly due to the in- 
tegrated effective path-length along the line of sight, rather 
than intrinsic variations. 

For the spectral line computation we use a gaussian line 
profile at a wavelength of 2.1218 fim with a cr = 
i.e. a FWHM of 33 km s“^, which ensures the profile is 
properly sampled with 128 channels over a full spectral range 
of ±700 km s~^ around the local standard of rest. 

Small-scale structure is added using various 
standard tools that modify the shape of the initially 
spherical 3D mesh. The general modeling proce- 
dure is described on the support website of SHAPE 
{www.astrosen.unam.mx/shape) . Modeling the loca- 
tion of individual small-scale features is based on 
the following assumption: the Homunculus can be 
approximated as a surface structure and that the im- 
age together with the expected velocity field yields a 
one-to-one mapping between the projected observed 
position and velocity to the position on the Ho- 
munculus surface (e.g. Sabbadin et al. 2004; Magnor 
et al. 2005, Steffen, Garcia-Segura Koning 2009) 

The model P-V diagrams (Fig. 3) and images have been 
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Figure 3. Observed and model P-V images. The observational P-V images (upper panel in each P-V image pair), selected from 
the approximately one hundred images available, are most representative of the spatial- velocity changes as the slit moves across the 
Homunculus along its main axis, with the slit direction perpendicular to the polar axis of the Homunculus. The offset position of the 
numbered slit from the central source, in arcseconds, is indicated for each pair of P-V images (see Fig. 1). The spatial and spectral ranges 
of the P-V diagrams are —6.93" to +6.93" and —750 km s“^ to +800 km s“^, respectively. 


individually normalized and use a linear grayscale coding. 
As this paper focuses on the structure of the Homunculus, 
the overall relative brightness distribution is used only as a 
general guideline. 


4 RESULTS 

In this section we describe the results of modeling the small- 
scale structures of the bipolar Homunculus using the H 2 
observations presented in this paper, as well as some basic 
results regarding the location of the scattering dust relative 
to the emitting H 2 . For simplicity, we use the terms ‘blue 
lobe’ and ‘red lobe’ when referring to the Homunculus’s ap- 
proaching (SE) and receding (NW) lobes, respectively. 


4.1 Velocity field 

As mentioned, the assumption of a homologous expansion 
for the Homunculus leads to a reconstruction of the struc- 
ture along the line of sight that is inconsistent with the con- 
tinuum image and the expectation that the molecular line 
emission originates outside of the dust. We find that within 
the precision of the data and the reconstruction method, a 
non-homologous expansion is only required for the blue lobe. 
For the blue lobe, the model includes a non-homologous ve- 
locity field in the form of a small poloidal velocity component 
with a latitudinal variation (Steffen, Garcfa-Segura & Kon- 
ing 2009). Here the velocity field is divided into two locally 
perpendicular components, i.e. a radial and poloidal velocity. 
The radial velocity Vr{r) = 680 km (r/3.4 x 10^^ m)^'^ 
is a function that is slightly non-linear, increasing faster 
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Figure 4. The 3D mesh structure of the SHAPE Homunculus model. For clarity, the mesh resolution has been reduced. Left: Orientation 
as viewed from Earth, with North up and East to the left. The directions, as projected onto the equatorial plane, of orbital apastron and 
the new near-equatorial kinematic features (e.g. the Red and Blue Protrusions) are marked with grey lines. Right: Similar to left, but 
for a non-inclined, non-rotated view where the polar axis is parallel to the plane of the page, showing that the Red and Blue Protrusions 
extend from the NW-red and SE-blue lobes, respectively, and do not lie in the equatorial plane. 


than linearly with distance r from the central star. The 
poloidal component V(p{(f>) has a broad peak at mid an- 
gles 45°) between the equator and pole of the blue 
lobe, going slowly to zero at the equator and pole. The 
peak poloidal velocity amounts to ~ 50 km s“^ and is 
hence only of order 10% of the peak expansion velocity 
of the Homunculus. For the red lobe the radial velocity is 
Vr{r) = 600 km s“^ (r/3.4 x 10^^ m) and hence linear, but 
slightly smaller than that of the blue lobe. No poloidal com- 
ponent was necessary for the red lobe. We hnd that these 
small corrections solve the problem of different reconstruc- 
tions based on images or Doppler-measurements, yielding a 
large-scale cylindrical outline that is consistent with both 
types of observations. 

The precision of the exponent 1.2 in the radial 
component and a peak ~ 50 km s~^ for the poloidal 
components of the velocity in the blue lobe is uncer- 
tain. These values should be taken only as an clear 
indication that there is some large-scale deviation 
from a homologous expansion, if the systematic de- 
viation between Doppler-velocity and image based 
reconstructions is not real. 

In the kinematic analysis of hydrodynamic simulations 
of wind-driven bipolar nebulae presented by Steffen, Garcfa- 
Segura & Koning (2009), the poloidal component points 
away from the pole, while here we hnd that the direction 
is toward the pole, implying a recohimation. While the ef- 
fect is only a small percentage of the total velocity, it is 


large enough to change the overall shape of the reconstructed 
Homunculus. While the study by Steffen, Garcia- Segura & 
Koning (2009) was aimed at planetary nebulae driven by a 
continuous stellar wind, in Eta Garinae an explosive expan- 
sion and different environmental conditions might produce 
the inverted poloidal velocity component. Such a how is ob- 
tained for instance in a Ganto-Flow, where a wind is shocked 
on the inside of a lobe in an external density gradient and is 
redirected toward the pole of the lobe (Ganto & Rodriguez 
1980). New hydrodynamical simulations might shed light on 
the mechanisms responsible for this deviation from homolo- 
gous expansion, but are outside the scope of this paper and 
must be left for future research. 


4.2 Small-scale structures 

The Homunculus is highly structured on small scales at 
visible and IR wavelengths, and departures from axisym- 
metry, especially in the blue lobe, have been known for 
some time (Morse et al. 1998; Smith 2006). Filamentary 
dust structures highlight the observed speckled structure 
of the lobes (Fig. 1). Furthermore, there are apparently 
radially-expanding flat spikes in the equatorial plane and 
polar “holes” in the lobes (Smith, Gehrz & Krautter 1998). 

The identihcation of symmetric small-scale structures in 
the Homunculus is expected to provide clues to the ejection 
mechanism in relation to the existence of a highly eccentric 
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Feature 

Longitude 

A Longitude 

Latitude 

A Latitude 

Blue trench 

180 

130 

75 

10 

Red trench 

0 

130 

75 

10 

Blue “hole” 

0 

N/A 

8 

15 

Red “hole” 

180 

N/A 

8 

12 

Blue Protrusion 

-55 

30 

10 

7 

Red Protrusion 

55 

30 

10-30 

20 


Table 1. Summary of the approximate geometric properties of the small-scale features identified in the 3D model. Angles are measured 
in degrees. Longitude is with respect to the apastron direction of the binary star and positive is counter-clockwise from Earth’s point of 
view. Here A refers to the angular spatial extent of a particular feature, not the uncertainty in the measurement. The uncertainties in 
the angular values are estimated to be of order 10 — 20%. In the case of the latitude of the polar holes, the latitudinal angle is given as 
the deviation from the main axis of the lobes in the model. The uncertainty of this value is ^ 2 ° . 



regularities 


Red Protrusion 




B me Protrusion 


Blue trench 


Blue lobe 


Blue Protrusion 


Red trench 


Figure 5. Shaded 3D view of the SHAPE model Homunculus 
with the view from Earth at the top (North is up. East is left) 
and a ‘flipped’ view from the opposite direction at the bottom 
(North is down. East is still left). Click image for a 3D interactive 
view (Adobe Reader® only). 


central binary with interacting winds. We thus attempt to 
model small-scale features in the observed H 2 P-V diagrams 
in some detail. The features that we find with blue and red 
counterparts are identified in Fig. 5 and their approximate 
geometric properties are summarized in Table 1. 

The small-scale features have been modeled in- 
teractively by visual superposition and comparison 
in SHAPE. While in principle possible (Nakashima 


et al., 2010), automated reconstruction of high- 
contrast features such as those found in our obser- 
vations will not work appropriately. The procedure 
is also limited to radial outflows. Due to the large 
amount of data and parameters to fit, a formal least- 
squares type fitting procedure is impractical. SHAPE 
can do this for final adjustments of inclination an- 
gles and similar parameters. But the result is usu- 
ally not significantly better than an interactive vi- 
sual approach. Tests with similar features placed at 
random positions on an artificial objects similar to 
the Homunculus and then recovering the position 
and structure revealed a precision of azimuthal po- 
sition of 5—10°. 

High contrast display of some of the P-V diagrams re- 
veals strong deviations from a smooth lobe structure near 
the equatorial region (the features labeled as the Blue and 
Red Protrusions in Fig. 3). Our modeling shows that these 
are low-brightness protrusions from the main lobes located 
at highly symmetric positions with respect to the orienta- 
tion of the binary orbit (left panel of Fig. 4). While the 
detailed structure of the Blue and Red Protrusions is differ- 
ent from each other, they straddle the direction of binary 
apastron by angles of approximately 55° on either side, as 
projected onto the orbital plane. The Blue and Red Protru- 
sions appear to be on the approaching and receding sides 
of orbital periastron, respectively (see e.g. fig. 12 of Madura 
et al. 2012). The modeling shows that the protrusions are 
not part of any equatorial outflow, they are clearly outside 
of this plane and part of the lobes (right panel of Fig. 4 and 
Fig. 5). We find that the Blue Protrusion is approximately 
10° out of the equatorial plane, while for the Red Protrusion 
a specific direction cannot be determined accurately as it is 
spread over a range of latitudes from ~ 10° to 30°. 

In the polar regions there are two clearly symmetric 
features, in addition to smaller “irregularities” of which the 
symmetry could not securely be determined. The symmet- 
ric features are the “polar holes” and what we call “polar 
trenches” (Fig. 5). The polar holes near the centres of the 
lobes are well known from spectroscopic observations (e.g. 
Smith, Gehrz & Krautter 1998; Teodoro et al. 2008). How- 
ever, for the first time, we identify the trenches on the blue 
and red lobes as such. In optical observations the blue trench 
is easily visible as a continuous curved dark line crossing the 
polar region of the blue lobe. The red counterpart can only 
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be seen in spectroscopic observations that reveal the far side 
of the red lobe. 

The modeling reveals that the polar holes may not be 
completely void of dust or H 2 . Instead of actual holes, these 
features may be indentations. Since the wall on the far side 
of the blue hole is at a large angle to the line of sight, the op- 
tical depth is lower than in other regions of the Homunculus, 
allowing emission from the inside to escape more easily. Un- 
fortunately, determining whether or not these features are 
actual holes is difficult and depends on detailed knowledge of 
the exact spatial diameter of the holes, the thickness of the 
polar lobes, and the inclination of the Homunculus nebula 
(Teodoro et al. 2008). While we model the holes as inden- 
tations in Fig. 5, within the precision of the X-shooter data 
and reconstruction method, models that assume a hole with 
an inward border of finite thickness and models assuming an 
indentation without a hole both reproduce the observed P-V 
diagrams equally well. Determining whether these features 
are holes in the true sense rather than indentations will re- 
quire more precise measurements of the hole diameters and 
thicknesses of the polar lobes. Detailed investigation of the 
structure of the Homunculus in [Fe ii] would also help clarify 
the nature of the polar holes. 

The newly-identified trenches in the polar regions of 
the lobes are particularly striking features. They are point- 
symmetric to the central stars. The modeling indicates they 
are not quite half-circles, centered on the polar holes and en- 
compassing an angle of ~ 120° to 140°, similar to the angu- 
lar distance between the Blue and Red Protrusions discussed 
above. We cannot exclude the possibility that the trenches 
are full circles centered on the holes. The inclined nature of 
the Homunculus and our line of sight into the trench may 
produce the appearance of a partial trench. Even if this is 
the case, the modeling and 3D kinematic observations show 
that the trenches are point symmetric, with one side of the 
trench clearly deeper than the other. 

As was the case with the polar holes, it is difficult 
to tell whether the trenches are complete rings without 
higher-resolution data. The available data and models can- 
not determine whether the trenches form complete circles. 
Whether there is any physical significance in these unexpect- 
edly highly symmetric features is unclear at this stage. The 
point-symmetric nature of the trenches implies they may 
be unrelated to the highly eccentric central binary and are 
possibly a result of the originating explosion. The red and 
blue trenches have the same angular extent and appear to 
be symmetric around the apastron direction. This may be 
a coincidence or it may indicate some connection between 
these features and the central binary. How these structures 
might be related to circumstellar matter existing prior to 
the event that led to the formation of the Homunculus, or 
the ejection event itself, is an intriguing question for future 
observational and theoretical work. 


4.3 Effects of H 2 and dust mixing 

In order to better understand the observed brightness dis- 
tribution in the H 2 data cube and its representation as P-V 
diagrams, we have added a dust component to the Homuncu- 
lus and computed the absorption and scattering of H 2 with 
radiation transfer along the line of sight. Before introducing 
inhomogeneities in the H 2 emissivity distribution, we first 




Figure 6. Rendered H 2 model images with (top) and without 
(bottom) dust. The top view assumes the H 2 and dust are fully 
mixed. Since the dust-to-emission ratio decreases with distance, 
the image without dust is noticeably brighter in the central region. 


investigated what brightness changes may be due to dust 
absorption. Therefore we used a smooth H 2 emissivity and 
dust distribution that decrease in density with distance r 
from the central star as oc (ro/(r + ro/5))“, with a = 2. This 
expression has been chosen as an approximation to 
geometric expansion with distance assuming and ini- 
tially uniform density, ro is a reference distance and 
ro/5 is used to avoid the expression to tend to infin- 
ity for small distances. 

The dust density was set high enough to partially ab- 
sorb the H 2 emission. We used a silicate dust radiation trans- 
fer model as described in Koning et al. (2013) with the dust 
properties listed in Section 3.1. Isotropic scattering is as- 
sumed. Fig. 6 shows the final rendered H 2 images with (top) 
and without (bottom) dust. In the top image, the dust and 
H 2 are assumed to be mixed. 

To verify the expectation that the H 2 emission is mixed 
with the protective dust, with a thin layer of H 2 outside the 
dust shell, we considered three different cases for the relative 
spatial distribution of dust and molecular hydrogen. Case 1 
is that of complete mixing of both components. Cases 2 and 
3 are those of the dust being located mostly inside or outside 
of the molecular shell, respectively. For Cases 2 and 3 we 
generated a copy of the H 2 mesh, assigned a dust species to 
it and scaled it in such a way that the dust is either inside or 
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Figure 7. Example SHAPE model P-V diagrams illustrating the 
three different cases investigated for the relative spatial distribu- 
tion of dust and H 2 in the Homunculus lobes. The chosen P-V 
diagram corresponds to slit position 30 in Eig. 3, which is for a 
slice through the middle of the blue lobe (Eig. 1). Top left: Dust 
is located primarily interior to the H 2 shell. Top right: Dust and 
H 2 are completely mixed. Bottom left: Dust is located primarily 
exterior to the H 2 shell. Bottom right: Observed P-V diagram. 


outside the H2 shell. The density of the dust was adjusted 
such that for each case the best match was found. 

A comparison of the three cases (Fig. 7) shows that 
the P-V diagrams for Cases 1 and 3 deviate considerably 
from the observations, while Case 2 is able to reproduce the 
most salient brightness variations. As expected and observed 
in the P-V diagrams, the brightness of the H2 emission on 
the backside of the P-Vs is diminished (see also Fig. 3). 
Some of the gaps in the P-V diagrams of the backside are 
reproduced and, in the model, are caused by increased dust 
absorption due to higher column densities. This is an effect 
of the local geometry, rather than intrinsic density variations 
of H2. Here the viewing direction is roughly tangential to the 
surface yielding strong absorption of the far side. Most P-V 
images therefore show a reduced brightness near the top and 
bottom edges of the velocity ellipses. 

When the dust is located outside of the molecular shell, 
both sides of the H2 Homunculus receive the same level of 
absorption from only the front side of the dust shell, and 
there is not much difference in brightness between the fronts 
and backs of the lobes (bottom left panel of Fig. 7). When 
the dust is located interior to the molecular shell, only the 
backside of the H2 shell is absorbed, by both the front and 
back of the dust shell. This produces the strongest contrast 
between the fronts and backs of the lobes in the P-V dia- 
grams (top left panel of Fig. 7). When the dust and H2 are 
fully mixed, the front is partially absorbed by only the front 
dust, while the back side is fully absorbed by the front dust 
and partially by the back dust, fn reality, the situation is 
likely between one where the dust is located interior to the 
H2 shell and fully mixed with it. 


Strong dust absorption leads to strong variations in the 
P-V brightness, which, on the scale of these observations, 
have been reproduced by the model without intrinsic small- 
scale variations of the emissivity. The P-V diagrams still 
have many additional brightness variations though, due to 
intrinsic local emission and absorption variations which have 
not been modeled. 


5 DISCUSSION 

The improved 3D model of the Homunculus presented above 
provides new insights on the detailed smaller-scale structure 
of the bipolar lobes. Two of the most significant new re- 
sults of this work are the identification of the protrusions 
located near the equatorial plane in both lobes and the 
spatially-extended polar trenches, fnterestingly, both fea- 
tures are symmetric about the apastron direction of the bi- 
nary orbit, with the angular distance between the Blue and 
Red Protrusions 110°) being very similar to the angu- 
lar extent of each polar trench 130°). We point out that 
these angles are nearly identical to the total opening angle 
in the orbital plane of the current wind-wind collision cav- 
ity, which is ~ 110° when using a secondary /primary wind 
momentum ratio 77 0.12 (see Madura et al. 2013). 

The similarity in value between the colliding- 
winds-cavity opening angle and the angular separa- 
tion/extent of the newly discovered small-scale fea- 
tures (i.e. the blue and red protrusions and polar 
trenches) may be a numerical coincidence. However, 
a key point to keep in mind regarding the Eta Cari- 
nae binary is the systems extremely high eccentric- 
ity 0.9), which causes the companion star to spend 
most of the 5. 54- year orbit on the apastron side of 
the system. This creates a large-scale (> 1600 au) 
time-averaged wind cavity on the apastron side of 
the system whose total opening angle is again ~ 110° 
at nearly all phases of the orbit, except for the short 
time near periastron passage when the companion 
passes behind the primary (see e.g. figures 8 and B5 
of Madura et al. 2013). Interestingly, for the binary 
orientation derived by Madura et al. 2012, this time- 
averaged cavity structure lies near the Homuncu- 
lus equatorial plane and along those directions on 
the sky where the red and blue protrusions are ob- 
served. The similarity of the time-averaged wind- 
wind-cavity opening angle and the angular extent of 
the polar trenches is more likely a coincidence given 
their apparent point-symmetric nature. 

When projected onto the orbital plane, the Red and 
Blue Protrusions appear to be located on the apastron side 
of the orbit. While the polar holes (Fig. 5) are not a new dis- 
covery, we find that both appear to be symmetrically offset 
from the polar axis of the large-scale Homunculus by ~ 8° 
and aligned with the projected apastron direction of the or- 
bit. However, we note that the orbital axis of the binary is 
also slightly offset from the polar axis of the Homunculus 
and nearly aligned with the projected (on the sky) direction 
of apastron ((Madura et al. 2012)). Therefore, the holes may 
be more directly related to the binary orbital axis than the 
projected direction of apastron, especially since the orbital 
axis and polar holes are closely aligned in 3D space. Alter- 
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natively, the polar holes may be aligned with the rotation 
axis of the central star, assuming that the stellar-rotation 
and binary orbital axes are also closely aligned. 

With the exception of the smaller polar “irregularities” 
(Fig. 5), which are likely a result of various instabilities in 
the thin Homunculus shell (see e.g. Smith 2013) and/or in- 
teraction of the expanding Homunculus with pre-existing 
material, all small-scale features mentioned appear to be 
aligned relative to the apastron direction of the central bi- 
nary (although in the case of the point-symmetric polar 
trenches this may be just a coincidence). We therefore con- 
clude that the interaction between the outflows and/or radi- 
ation from the central binary stars and their orientation in 
space has had, and possibly still has, a strong influence on 
the structure of the Homunculus nebula. Such interactions 
might help explain the non-homologous expansion and small 
poloidal velocity component required to properly model the 
blue lobe (Section 4.1). 

Of particular relevance is the fast, lower-density wind 
from the companion star, which at Uoo ~ 3000 km s“^, is 
~ 5 to 40 times faster than the expanding Homunculus, de- 
pending on the latitude. In the ~ 170 years since the Great 
Eruption, the fast secondary wind has had enough time to 
catch up to and interact with portions of the Homunculus 
expanding in our direction, most especially at low latitudes 
near the equator where the expansion of the Homunculus 
is slowest (Smith 2006). Portions of the Homunculus fac- 
ing away from ns are mainly unaffected by the companion 
star’s wind since the much denser wind of the LBV primary 
is in the way and the colliding wind cavity opens in our 
direction. We speculate that some of the observed symmet- 
ric small-scale features in the Homunculus (more speciflcally 
the Red and Blue Protrusions) may be a result of interac- 
tions with the fast secondary wind. Detailed 3D numerical- 
hydrodynamical simulations beyond the scope of this paper 
are necessary to test this hypothesis, but such an interaction 
may help explain another peculiar feature of the Homuncu- 
lus, the thin equatorial skirt. 

Just as Smith (2006) noted, we And no detectable NIR 
H 2 emission from the thin equatorial skirt that is promi- 
nent in optical images of the Homunculus, implying that the 
skirt is of much lower density than the walls of the bipolar 
lobes and contains signiflcantly less mass, < 0.5 Mq (Smith 
2006). The extremely fast and low-density secondary wind 
has likely collided with the near-equatorial regions of the 
expanding Homunculus. The large difference in velocity be- 
tween the secondary wind and Homunculus should result in 
various instabilities at the collision interface (e.g. Rayleigh- 
Taylor, Kelvin-Helmholtz, non-linear thin shell). This, com- 
bined with the lack of any dense material exterior to the Ho- 
munculus walls, should disrupt the thin Homunculus shell 
near the equatorial region, allowing the wind of the sec- 
ondary to eventually plough through it. 

The above-described situation is analogous to that ob- 
served in 3D smoothed particle hydrodynamics (SPH) simu- 
lations of Eta Carinae (see e.g. flg. 8 of Madura et al. 2013). 
These simulations show that while the companion and its 
wind are at periastron, a dense, thin shell of primary wind 
flows toward the observer. When the binary companion re- 
turns to the apastron side, its fast wind collides with this 
dense shell of primary wind. Eventually, within half of an or- 
bital cycle, instabilities in the shell cause it to break apart. 


allowing the fast secondary wind to plough through it. This 
results in a ‘spray’ of blobs of primary wind moving in a 
sea of fast, low-density secondary wind, concentrated in the 
orbital plane. The skirt of the Homunculus may likewise be 
the remains of the fast secondary wind ploughing through 
the thin dense shell of the Homunculus, forming a spray of 
material concentrated in the equatorial region. 

If the secondary’s wind can escape from the Homuncu- 
lus’s equatorial regions, then so should some of the sec- 
ondary’s ionizing radiation. However, the secondary’s ra- 
diation would not escape uniformly in all directions, but 
would escape preferentially in directions of lowest density, 
with some regions/directions shielded from the hardest UV 
flux by the remaining denser blobs of material that previ- 
ously composed the walls of the Homunculus. The combi- 
nation of ionizing secondary radiation and low-density sec- 
ondary wind in the skirt could explain why no NIR H 2 emis- 
sion is observed there. Moreover, because this process occurs 
only on the apastron side of the binary, the equatorial skirt 
would appear asymmetric and localized to the side of the 
Homunculus facing ns, just as observed in optical images. 

We add that this interpretation for the formation of the 
equatorial skirt could help resolve another of its unexplained 
peculiarities, namely, the presence of younger ejecta from the 
1890 eruption that appears to coexist with older ejecta in 
the skirt (Smith & Gehrz 1998; Davidson et al. 2001; Smith 
et al. 2004; Smith & Townsend 2007). In order to continue 
ploughing through the outer, more massive Homunculus, the 
fast secondary wind would need to penetrate the inner Lit- 
tle Homunculus, which is thought to have formed during the 
1890 eruption (Ishibashi et al. 2003; Smith et al. 2004). Over 
time, the fast secondary wind and remnants of the equato- 
rial regions of the Little Homunculus shell would be driven 
outward. Eventually, the blobs of older skirt material and 
newer material from the Little Homunculus would mix with 
the secondary wind. We therefore suggest that the equato- 
rial skirt, and possibly the Red and Blue Protrusions, may 
be a result of the interaction of the fast secondary wind with 
the Homunculus and Little Homunculus post-eruption. 

Of course another possibility is that the symmetric 
small-scale features in the Homunculus are intrinsic to the 
Great Eruption. The strongest evidence for this possibility 
comes from the point-symmetric nature of the polar trenches 
and polar holes. Even if this is the case, the relation of these 
symmetric features to the apastron direction of the orbit 
and/or the binary orbital axis strongly implies that bina- 
rity had some role in helping trigger the eruption and/or 
shape the nebula. More evidence for a link between bina- 
rity and the Great Eruption comes from the correlation of 
periastron passage with the brightening peaks observed in 
Eta Garinae’s light curve in the decades leading up to the 
eruption (see flg. 2 of Damineli 1996 and also Smith & Erew 
2011). A stellar collision scenario such as that envisioned 
by Smith (2011) or a binary interaction like that in Kashi 
& Soker (2010) may have occurred during the eruption, al- 
though the details are far from certain. 

Another pair of peculiarities worth discussing are the 
polar holes of the Homunculus. Are these true holes or just 
indentations in the lobes? Were the holes created during 
the original eruption event, or signiflcantly later afterward? 
What physical process formed them? The holes may be the 
result of instabilities in the expanding bipolar structure, sim- 
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ilar to those that likely created the smaller irregularities seen 
in the NW red lobe. However, one would have to explain why 
such instabilities produced a pair of point-symmetric holes 
that are offset from the polar symmetry axis of the neb- 
ula. Instead, perhaps some sort of polar jet closely aligned 
with either the central star’s rotation axis or the binary’s 
orbital axis was present and produced the holes. This jet 
may have been present during the initial eruption or created 
much later. Unfortunately, we do not have enough informa- 
tion to constrain the timing of any hypothetical jets, nor can 
we currently explain the potential mechanisms behind their 
formation (see e.g. Soker 2003). Detailed examination of the 
Homunculus in [Fe ii] emission and of the interior Little Ho- 
munculus, searching for the presence of similar polar holes, 
would help greatly in understanding their origin. 

One final note regarding the polar holes. It would be 
interesting to investigate what the existence of these fea- 
tures imply for theories for the formation of the Homunculus 
that rely on mass-loss from a rapidly-rotating star (Owocki 
& Gayley 1997; Owocki et al. 1998; Dwarkadas & Owocki 
2002; Owocki 2003; Smith 2006). In such models, the mass- 
loss, and hence density, is expected to be highest over the 
poles. Yet, the mass fraction as a function of latitude in the 
lobes clearly shows a peak between 50° and 60° (see fig. 5b 
of Smith 2006), with a much smaller mass fraction at the 
poles. Combined with the difficulties of a steady radiation- 
driven wind powering such extreme mass-loss (Owocki et 
al. 2004; Smith & Townsend 2007), this may indicate that 
steady mass-loss from a rapidly rotating single star is not 
the correct explanation. 

The newly-discovered point-symmetric polar trenches 
and near-equatorial Red and Blue Protrusions, as well as the 
polar holes, pose some interesting complications for models 
for the Homunculus that are based on an explosion from a 
rapidly-rotating single star. Prior to the observations and 
3D model in this paper, the model by Smith & Townsend 
(2007) seemed a potentially viable explanation for the shape 
of the Homunculus. However, that model produces an ax- 
isymmetric nebula and skirt that do not match the more 
detailed structures apparent in the new observations pre- 
sented above. We emphasize that this does not rule out the 
possibility that the central primary star was near critical ro- 
tation for mass loss or that some sort of explosion occurred 
which helped power the Great Eruption. Instead, it more 
strongly implies that single-star models for the formation 
and shape of the Homunculus are not correct, and that bi- 
narity and interacting stellar winds played an important role 
in shaping the detailed 3D morphology of the Homunculus. 
It would be interesting to see if the inclusion of binarity 
and colliding winds in a model similar to that of Smith & 
Townsend (2007) can produce the more detailed 3D struc- 
ture of the Homunculus. 

As one can clearly see, the Homunculus is an ex- 
tremely complex structure whose formation and shape will 
likely only be adequately explained with the assistance of 
3D numerical-hydrodynamical simulations. Numerous pos- 
sible scenarios exist for triggering the eruption and shaping 
the large- and small-scale features observed. Future models 
should strongly consider the effects of binarity, collisions, 
and the spatially-extended colliding stellar winds. Stellar 
mergers should also be explored. The pre-emption envi- 
ronment immediately surrounding the central stars is also 


expected to be of some importance in influencing the re- 
sultant asymmetric shape and various possible instabilities, 
and should be included if possible. 


6 CONCLUSIONS 

In this paper we presented new kinematic molecular hydro- 
gen data from the ESO VLT/X-shooter instrument that 
fully map the Homunculus nebula around Eta Garinae. 
Based on these observations we constructed the first full 
3D model of the shape of the bipolar nebula that includes 
small-scale structures. As in Smith (2006), we find no de- 
tectable NIR H 2 emission from the thin equatorial skirt that 
is prominent in optical images of the Homunculus. By intro- 
ducing a non-radial correction to the velocity field, we find 
a unique solution for the large-scale structure of the bipolar 
Homunculus, unifying earlier attempts that yielded different 
solutions for models based on imaging and kinematic data. 

We further extended the model to incorporate observed 
small-scale deviations from axisymmetry. The geometric po- 
sitioning of these features is summarized in Table 1. Based 
on our observations, we accurately place the positions of the 
well-known polar holes, finding that they appear to be sym- 
metrically offset from the axis of the large-scale Homunculus 
by approximately 8°. In both cases this deviation is closely 
aligned with the apastron direction of the central binary or- 
bit and the 3D orientation of the binary’s orbital axis. 

Similarly, we find that the newly identified polar 
trenches (Fig. 5) are point-symmetric with respect to the 
central stars and may also be related to the apastron direc- 
tion of the binary orbit. The most intriguing new features 
are the Blue and Red Protrusions in the lobes (Figs. 3-5). 
They are not exactly in the equatorial plane, but of order 
10° in latitude out of the plane. In longitude they are both 
about 55° away from the direction of apastron, in opposite 
directions. It might be noteworthy that the angular distance 
between the Blue and Red Protrusions (~ 110°) is similar 
to the angular extent of each of the polar trenches 130°), 
and nearly identical to the estimated total opening angle of 
the current-day inner wind- wind collision cavity (~ 110°). 

The smaller “irregular” features in Fig. 5 have no ap- 
parent symmetry or relation to the orientation of the bi- 
nary orbit. They might be attributable to instabilities in 
the expanding Homunculus. All other major small-scale fea- 
ture appear to be aligned relative to the apastron direc- 
tion of the central stars, although in the case of the point- 
symmetric polar trenches this may be a coincidence. We 
therefore conclude that the interaction between the outflows 
and/or radiation from the central binary stars and their ori- 
entation in space has had, and possibly still has, a strong 
influence on the structure of the Homunculus nebula. In- 
teraction of the fast wind from the secondary star with the 
near-equatorial regions of the Homunculus and Little Ho- 
munculus may help explain the protrusions and thin equa- 
torial skirt. The possibility that the current stellar winds 
are influencing the Homunculus’s structure, plus the require- 
ment of a non-homologous expansion to properly model the 
SE blue lobe, imply that the standard assumption of a ho- 
mologous expansion may no longer hold. 

Our new observations and 3D model pose serious com- 
plications for single- star theoretical models of the Homnncu- 
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lus’s formation. Prior to this work, there was no strong rea- 
son to consider a binary model for the formation or shape 
of the Homunculus. Detailed 3D numerical- hydrodynamical 
modeling of the Homunculus’s formation that includes the 
central binary and its interacting winds are now required to 
determine if such interactions can simultaneously reproduce 
the newly-identified small-scale features, overall large-scale 
bipolar shape, the polar holes, and thin equatorial skirt. 


ACKNOWLEDGMENTS 

W. S. acknowledges financial support through grant UNAM- 
PAPHT IN101014. M. T. is supported by CNPq/MCT- 
Brazil through grant 201978/2012-1. T. I. M. is supported 
by an appointment to the NASA Postdoctoral Program 
at the Goddard Space Flight Center, administered by 
Oak Ridge Associated Universities through a contract with 
NASA. J. H. G. is supported by an Ambizione fellowship 
of the Swiss National Science Foundation. A. D. acknowl- 
edge FAPESP for continuous financial support. The authors 
thank D. Clark for useful discussions. 


REFERENCES 

Allen, D., & Hillier, J. 1991, Proceedings of the Astronom- 
ical Society of Australia, 9, 120 
Allen, D. A., & Hillier, D. J. 1993, Proceedings of the As- 
tronomical Society of Australia, 10, 338 
Canto & Rodriguez, 1980, ApJ, 239, 982-987 
Chesneau, O., Min, M., Herbst, T., et al. 2005, A&A, 435, 
1043 

Clark, D. M., Lopez, J. A., Steffen, W., & Richer, M. G. 
2013, AJ, 145, 57 

Cox, P., Mezger, P. G., Sievers, A., Najarro, F., Bronfman, 
L., Kreysa, E., & Haslam, G. 1995, A&A, 297, 168 
Currie, D. G., Dowling, D. M., Shaya, E. J., et al. 1996, 
AJ, 112, 1115 

Damineli, A. 1996, ApJ, 460, L49 

Damineli, A., Conti, P. S., & Lopes, D. E. 1997, New As- 
tronomy, 2, 107 

Davidson, K. & Humphreys, R.M. 1997, ARA&A, 35, 1 
Davidson, K., & Humphreys, R. M. 2012, Astrophysics and 
Space Science Library, Vol. 384, Eta Carinae and the Su- 
pernova Imposters 

Davidson, K., Smith, N., Gull, T. R., Ishibashi, K., & 
Hillier, D. J. 2001, AJ, 121, 1569 
Dwarkadas, V. V., & Balick, B. 1998, AJ, 116, 829 
Dwarkadas, V. and Owocki, S.P. 2002, ApJ, 581, 1337 
Erank, A., Balick, B., & Davidson, K. 1995, ApJ, 441, L77 
Erank, A., Ryu, D., & Davidson, K. 1998, ApJ, 500, 291 
Gallagher, J. S. 1989, lAU Colloq. 113: Physics of Lumi- 
nous Blue Variables, 157, 185 
Garcia-Diaz, M.T., Lopez, J.A., Steffen, W., Richer, M.G., 
2012, ApJ, 761, 172 
Gaviola, E. 1950, ApJ, 111, 408 

Goldoni, P., Royer, E., Erangois, P., et al. 2006, Ground- 
based and Airborne Instrumentation for Astronomy. 
Edited by McLean, 6269, 80 

Gomez, H. L., Vlahakis, C., Stretch, C. M., et al. 2010, 
MNRAS, 401, L48 


Gonzalez, R. E., de Gouveia Dal Pino, E. M., Raga, A. C., 
& Velazquez, P. E. 2004a, ApJ, 600, L59 
Gonzalez, R. E., de Gouveia Dal Pino, E. M., Raga, A. C., 
& Velazquez, P. E. 2004b, ApJ, 616, 976 
Groh, J. H., Hillier, D. J., Madura, T. L, & Weigelt, G. 
2012a, MNRAS, 423, 1623 

Groh, J. H., Madura, T. L, Hillier, D. J., Kruip, C. J. H., 
& Weigelt, G. 2012b, ApJL, 759, L2 
Groh, J. H., Madura, T. L, Owocki, S. P., Hillier, D. J., & 
Weigelt, G. 2010a, ApJL, 716, L223 

Groh, J. H., Nielsen, K. E., Damineli, A., et al. 2010b, 
A&A, 517, A9 

Hackwell, J. A., Gehrz, R. D., & Grasdalen, G. L. 1986, 
ApJ, 311, 380 

Hillier, D. J., & Allen, D. A. 1992, A&A, 262, 153 
Hillier, D. J. 1997, Luminous Blue Variables: Massive Stars 
in Transition, 120, 287 

Hillier, D. J., Davidson, K., Ishibashi, K., & Gull, T. 2001, 
ApJ, 553, 837 

Hillier, D. J., Gull, T., Nielsen, K., et al. 2006, ApJ, 642, 
1098 

Humphreys, R. M. & Davidson, K. 1994, PASP, 106, 1025 
Iben, L, Jr. 1999, Eta Carinae at The Millennium, 179, 367 
Ishibashi, K., et al. 2003, AJ, 125, 3222 
Kashi, A., & Soker, N. 2010, ApJ, 723, 602 
Koning, N., Kwok, S., & Steffen, W., 2013, ApJ, 765, 92 
Krugel, E., 2008, “An Introduction to Instellar Dust” , Pub- 
lisher: Taylor & Erancis, ISBN 9781584887072. 

Langer, N., Hamann, W.-R., Lennon, M., et al. 1994, A&A, 
290, 819 

Langer, N., Garcia-Segura, G., & Mac Low, M.-M. 1999, 
ApJ, 520, L49 

Lopez, J.A., Garcia-Diaz, M.T., Steffen, W., Riesgo, H., 
Richer, M.G., 2012, ApJ, 750, 131 
Madura, T. I. & Groh, J. H. 2012, ApJL, 746, L18 
Madura, T. L, Gull, T. R., Okazaki, A. T., et al. 2013, 
MNRAS, 436, 3820 

Madura, T. L, Gull, T. R., Owocki, S. P., et al. 2012, MN- 
RAS, 420, 2064 

Maeder, A., & Desjacques, V. 2001, A&A, 372, L9 

Magnor, M., Kindlmann, G., Hansen, C., Duric, C., 
2005, ’’Reconstruction and Visualization of Plan- 
etary Nebulae,” IEEE Trans. Visualization and 
Computer Graphics, vol. 11, no. 5, pp. 485-496 

Martin, J. C., Davidson, K., & Koppelman, M. D. 2006, 
AJ, 132, 2717 

Matt, S., & Balick, B. 2004, ApJ, 615, 921 
Meaburn, J., Wolstencroft, R. D., & Walsh, J. R. 1987, 
A&A, 181, 333 

Meaburn, J., Walsh, J. R., & Wolstencroft, R. D. 1993, 
A&A, 268, 283 

Meaburn, J. 1999, Eta Carinae at The Millennium, 179, 
159 

Mehner, A., Davidson, K., Eerland, G. J., & Humphreys, 
R. M. 2010, ApJ, 710, 729 

Morris, T., & Podsiadlowski, P. 2006, MNRAS, 365, 2 
Morse, J. A., Davidson, K., Bally, J., et al. 1998, AJ, 116, 
2443 

Morse, J. A., Kellogg, J. R., Bally, J., et al. 2001, ApJ, 548, 
L207 

Nakashima, J., Kwok, S., Zhang, Y., Koning, N., 2010, AJ, 
140, 490-499 


© 0000 RAS, MNRAS 000, 1-14 


14 W. Steffen et al. 


Okazaki, A. T., Owocki, S. P., Russell, C. M. P., & Corco- 
ran, M. F. 2008, MNRAS, 388, L39 
Ouyed, R., Kostka, M., Koning, N., Leahy, D.A., Steffen, 
W., 2012, MNRAS, 423, 1652 

Owocki, S. P., & Cayley, K. C. 1997, Luminous Blue Vari- 
ables: Massive Stars in Transition, 120, 121 
Owocki, S. P., Cayley, K. C., & Cranmer, S. R. 1998, Prop- 
erties of Hot Luminous Stars, 131, 237 
Owocki, S. P. 2003, A Massive Star Odyssey: From Main 
Sequence to Supernova, 212, 281 
Owocki, S. P., Cayley, K. C., & Shaviv, N. J. 2004, ApJ, 
616, 525 

Owocki, S. 2005, The Fate of the Most Massive Stars, 332, 
169 

Parkin, E. R., Pittard, J. M., Corcoran, M. F., & Ham- 
aguchi, K. 2011, ApJ, 726, 105 
Parkin, E. R., Pittard, J. M., Corcoran, M. F., Hamaguchi, 
K., & Stevens, 1. R. 2009, MNRAS, 394, 1758 
Pittard, J. M. & Corcoran, M. F. 2002, A&A, 383, 636 
Polomski, E. F., Telesco, C. M., Pina, R. K., & Fisher, R. S. 
1999, AJ, 118, 2369 

Ribeiro, V.A.R.M., Bode, M.F., Darnley, M.J., Barnsley, 
R.M., Munari, U., Harman, D.J., 2013, MNRAS, 433, 
1991-1996 

Sabbadin, F., Turatto, M., Cappellaro, E., Benetti, S., 
Ragazzoni, R., 2004, A&A, 416, 955-981 
Smith, C. H., Aitken, D. K., Moore, T. J. T., et ah 1995, 
MNRAS, 273, 354 
Smith, N. 2002, MNRAS, 337, 1252 

Smith, N., Cehrz, R. D., Hinz, P. M., et ah 2003a, AJ, 125, 
1458 

Smith, N., Davidson, K., Cull, T. R., Ishibashi, K., & 
Hillier, D. J. 2003b, ApJ, 586, 432 
Smith, N. 2006, ApJ, 644, 1151 
Smith, N. 2011, MNRAS, 415, 2020 
Smith, N. 2013, MNRAS, 429, 2366 
Smith, N., & Erew, D. J. 2011, MNRAS, 415, 2009 
Smith, N. & Townsend, R. 2007, ApJ, 666, 967 
Smith, N. & Davidson, K. 2001, ApJ, 551, LlOl 
Smith, N., & Cehrz, R. D. 1998, AJ, 116, 823 
Smith, N., Cehrz, R. D., & Krautter, J., 1998, AJ, 116, 
1332 

Smith, N. et ah 2004, ApJ, 605, 405 

Smith, N. & Owocki, S. R 2006, ApJ, 645, L45 

Soker, N. 2001, MNRAS, 325, 584 

Soker, N. 2003, ApJ, 597, 513 

Soker, N. 2004, ApJ, 612, 1060 

Soker, N. 2007, ApJ, 661, 490 

Steffen, W., Koning, N., Wenger, S., Morisset, C., Magnor, 
M., 2011, IEEE Transactions on Visualization & Com- 
puter Craphics, voh 17, num. 4, 454-465 (The arXive ref- 
erence is 2010arXivl003.2012S). 

Steffen, W., Carcfa-Segura, C., Koning, N., 2009, ApJ, 691, 
696-704 

Thackeray, A. D. 1951, The Observatory, 71, 167 
Thackeray, A. D. 1956a, The Observatory, 76, 103 
Thackeray, A. D. 1956b, The Observatory, 76, 154 
Thackeray, A. D. 1961, The Observatory, 81, 99 
Teodoro, M., Damineli, A., Sharp, R.C., Croh, J.H., Bar- 
bosa, C.L., 2008, MNRAS, 387, 564-576 
Vernet, J., Dekker, H., Dodorico, S., et ah 2011, A&A, 536, 
105 


This paper has been typeset from a TgX / DTgX file prepared 
by the author. 


© 0000 RAS, MNRAS 000, 1-14 


